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ABSTRACT 

We have obtained long-slit, infrared spectra of Jupiter with the Anglo Australian 
Telescope in the K and H bands at a resolving power of 2260. Using a line-by-line, 
radiative transfer model with the latest, improved spectral line data for methane and 
ammonia, we derive a model of the zonal characteristics in the atmosphere of this giant 
planet. We fit our model to the spectra of the zones and belts visible at 2.1 /mi using 
different distributions of cloud opacities. The modeled spectra for each region match 
observations remarkably well at K band and in low pressure regions at the H band. 
Our results for the upper deck cloud distribution are consistent with previous models 
(Banfield et al.1998) fitted to low resolution, grism spectra. The ability to obtain 
and model high resolution planetary spectra in order to search for weakly absorbing 
atmospheric constituents can provide better constraints on the chemical composition 
of planetary atmospheres. 

Key words: Jupiter, radiative transfer, atmospheric effects, techniques: spectro- 
scopic, infrared 
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1 INTRODUCTION 

The infrared spectrum of Jupiter is dominated by absorption 
features due to methane (CH4) as well as other gases such 
as ammonia (NH3) and collision induced absorption due to 
H2 and He. Methane, however, has a very complex rovibra- 
tional spectrum and the available spectrum line data are 
limited. Many models for Jupiter's spectrum have therefore 
been based on low resolution absorption coefficients or k- 
distributi on data, rath e r than on li ne-by-line te c hniqu es. For 
example IStrong et ail l| 19931 ) and llrwin et alJ (2005) have 
presented k-distribution mod els for methane ab sorption at 
a resolution of 10 cm" , and Irwin et al l |l999h have pre- 
sented similar data for ammonia. Using data such as this 
it is possible to model spectra at relatively low resolution 
such as those from Galileo-NIMS. However, many infrared 
spectrographs on ground-based telescopes can provide much 
higher spectral resolution. Recently there have been signif- 
icant improvements in t he available spect r al line data for 
methane and ammonia (Alb ert et al . 2009: Rothman et al 



2010) 
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In this paper we present a set of observations of the 
spectrum of Jupiter at a spectral resolution of ~2260 and 
investigate our ability to model this data using line-by- 
line methods. We use the radiative transfer model VS- 
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TAR (Versatile Software for Transfer of Atmospheric Radia- 
tion) , which is specifically designed for this purpose l|Bailevl 
2006). It derives radiative transfer solutions for specified 
wavelength ranges in the transmission, emission or reflec- 
tion spectrum of a planet. This software is designed to 
be extremely versatile in the range of atmospheres it can 
be applied to. In the past we have used it for the atmo- 
spheres of Earth and Venus jBailev. Simpson fc Crisd 120071 : 
iBailev et aI1l200Sl : lBailevll2009f ) and recent additions to the 
system allow it to be used for the atmospheres of solar- 
system planets, exoplanets, brown dwarfs and cool stars. 

Methane is an important absorber, not just in Jupiter 
and the other solar system giant planets, but also in brown- 
dwarfs (the presence of methane absorption is the defining 
characteristic of a T-dwarf) and is likely to be important in 
many exoplanets. It h as been reported in the sp ectra of the 
transiting exoplanets jSwain et all l200si . l2009h . While hot 
methane presents additional issues compare d with the cool 
meth ane spectra required for Jupiter (e.g. iBorvsov et al.l 
2002) the problem of accurately modelling the methane ab- 
sorption spectrum is common to a full understanding of all 
these objects. Thus Jupiter can be considered a test object 
for our future ability to model the spectra of extrasolar plan- 
ets, as well as an interesting object in its own right. 

Despite decades of research dedicated to understand- 
ing the structure and composition of Jupiter's atmosphere, 
unanswered questions remain about the detailed distribu- 
tion and transport of gases over the surface of the most mas- 
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sive planet of the Solar system. The wealth of data across its 
optic al and infrared spectrum from ground-bas ed telescopes 
(e.g. IClark fc McCordl 1 19791 : lKarkoschkalll994h. Pioneer 1 
llTomasko et al.lll97sfT Voyag er (IConrath fe Gieraschlll986(h 



Galileo (ICarlson et alj 19961) and the C assini spacecraft 
i|Kunde et all 12004 iBellucci et all I2004T ) provided direct 
measurements of physical properties at different heights of 
the atmosphere, which can be used to constrain input pa- 
rameters in models. Spectral information is used not only 
to derive the atomic and molecular abundances of the at- 
mosphere given the physical properties such as temperature 
and pressure at different heights, but also to study the distri- 
bution of aerosols and clouds, and dynamic characteristics, 
which include wind and temperature patterns due to global 
weather effects. 

Jupiter is covered with permanent layers of clouds, 
which are arranged into longitudinal, relatively stable re- 
gions determined by atmospheric circulation. The bright and 
dark bands are called zones and belts respectively, which are 
understood to correspond to up- and down-welling of gas in 
the convective flow of the atmosphere. The most abundant 
constituents of Jupiter's atmosphere are hydrogen, helium, 
methane and ammonia. Direct measurements of the scatter- 
ing properties in Jupi ter atmosphere take n with the Galileo 
Probe Nephelometer l|Ragent et al.l fl99"H ) suggest that the 
dense clouds of condensates, mostly ammonia, form in the 
almost adiabatic troposphere above 1600 to 700 mb. In addi- 
tion the thinner, stratospheric hazes, which contribute to the 
overall atmospheric opacity, have been detected at heights, 
which vary depending on the longitude of the planet. For 
example they appear elevated in the polar regions, where 
hydrocarbons, such as C2H6 and C2H2 are suggested to 
exist as a product of photoche mical and auroral discharge 
processes (|Prvor fc Hordlll99ll) . Another molecule consid- 
ered in a compositio n of the strat osphere is photochemi- 
cally produced N 2 H 4 l|Strobellll98t) . The detailed composi- 
tion of stratospheric haze remains uncertain due to limita- 
tions of currently available data. A compre hensiv e revie w of 
Jupiter's cloud structure was presented by llrwinl (|l999T ). 

The observations described in this paper are long slit 
spectra at relatively high resolving power (R = A/ A A ~ 
2260) in the H and K band (Table [TJ. The data were ob- 
tained with the Infr ared Imager and Spectrograph 2 (IRIS2) , 
l|Tinnev et al.ll2004T ) at the Anglo- Australian 3.9 Telescope. 
The spectra were subdivided into latitudinal regions, that 
closely resemble the pattern of zones and belts visible in the 
acquisition image observed with CH4 filter at 1.673 /jm (Fig- 
ure [TJ. The details of calibration technique are discussed in 
Section [5] We used the VSTAR package to develop models 
of the cloud structure for each region to best fit the ob- 
served spectra. In Section [3] we describe the VSTAR model- 
ing strategy to achieve this goal, which includes the model 
parameters and the lists of spectral lines used for each con- 
stituent of the atmosphere. The derived cloud opacities at 
different latitudinal locations are presented in Section 3] and 
followed by discussion in Section [5] 



2 OBSERVATIONS AND DATA REDUCTION 

The images and the IRIS2 spectra of Jupiter were acquired 
on the 25th of May 2007 (Table [TJ. The long slit (7.7') was 




Figure 1. The narrow- band, acquisition image of Jupiter taken 
in the CH4 absorption region at 1.673 fim shows a characteristic 
band structure of clouds in the atmosphere of the planet. The 
black lines mark the borders between regions, which correspond 
to zones and belts defined by the level of absorption visible at 
2.1 lira through a 1" slit positioned vertically at the centre of the 
planetary disk. All observed spectra are derived as the averages 
over the number of pixels comprising each zone or belt. 



positioned across the meridian of the planet. The width of 
the slit is 1", which is equivalent to 2.2 pixels on the detector. 
The spectra cover the H (1.62-1.8 (im) and K (2.0-2.4 fj,m) 
bands. The narrow-band image of the whole planet in CH4 
(Figure [TJ shows the characteristic band cloud structure. 
The calibrated spectra were extracted and averaged over 
11 regions denoted by names given in Table [2] The region 
selection was based on the reflectivity profile averaged across 
the slit over the section of the K band spectrum in the region 
of the collision induced H2-H2 absorption between 2.08 and 
2.13 [im. The black horizontal lines overploted on the image 
in Figure [l] show the extent of each region. 

Data reduction was carri ed out using standard proce- 
dures of the Figaro package (|Shortridge et alH l995:). The 
wavelength calibration was achieved by using a Xenon lamp. 
The spectra were divided by spectroscopic flat fields in or- 
der to remove the non-uniform response of the detector. The 
normalized flat-field images were formed by subtracting the 
image fully illuminated by a quartz calibration lamp and the 
dark image taken with the lamp switched off. Sky subtrac- 
tion was achieved by taking four exposures with two spectra 
spatially offset in the same frame and subtracting two cor- 
responding frames with offset spectra from each other. 

Spectral curvature was removed by resampling the 
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frames using correction derived from spectra of standard 
stars taken for each band. Bad pixels due to cosmic rays 
and known imperfections of the detector were corrected by 
interpolation between their neighbours. 

We used standard stars, which were close in angular 
separation to the planet. They were chosen to match the 
solar spectral type as closely as possible. The Jupiter spec- 
trum was divided by the standard star spectrum to remove 
both telluric absorptions and solar features. The resulting 
spectrum corresponds to the r eflectance from the planet. 
iBailev. Simpson fc Crisd |2007l ) have pointed out potential 
problems with using standard stars for removal of telluric 
absorption for planetary spectra. However, in the case of 
Jupiter the absorption features in the planet are due to dif- 
ferent molecules from those in the terrestrial atmosphere 
spectrum and these problems are largely avoided. 

A true level of intensity can be derived by converting 
a number of counts, N c in the spectra into the flux density 
derived from the calibrated flux density, F st , of the standard 
star. We derive the radiance factor, l/Fsun, where I is the 
reflected intensity of the planet and irFsun is the solar flux 
incident on the planet from: 



Fsuri 



_ ±D 2 (N} upUer )F st 
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(1) 



The size of the slit, A s ut is given in radians, and the 
distance of Jupiter, D from the Sun is measured in astro- 
nomical units. The flux in each region is then scaled by the 
number of pixels, it encompasses. 

The typical spectrum and our VSTAR model in the 
equatorial region of Jupiter are presented in Figure [2] in 
the wavelength range between 1.62 and 2.4 microns. As 
expected the methane absorption clearly dominates in this 
spectral range. The H2-H2, collision induced smooth absorp- 
tion peaks at 2.12 |im, which appears to be very shallow in 
the polar region spectra. This reflects the contributions of in- 
creasingly slant paths through the atmosphere, which probe 
vertically higher regions with lower density of Eh. 

The moon, Europa was transiting in front of the planet 
during the observations and its reflected light was present in 
the slit while the spectra were taken. Its effect on the spectra 
of the Southern Equatorial Zone (SEZ) and Southern Belt 
(SB) is clearly be visible in Figure [7] showing the increased 
reflectivity in the regions of the largest methane absorption. 
Europa is thought to have the lar ge amount of pure wa ter ice 
on its surface (see the spectra in Calvin et al.l |l995)). The 



highly reflective shoulders of water bands observed in the 
spectra around 1.8 and 2.2 /jm contribute to the combined 
high albedo of the planet and the moon in the regions where 
otherwise the planetary methane absorption is the highest. 



3 THE VSTAR MODELING PROCEDURE 

The Versatile Software for Transfer of Atmospheric Radia- 
tion (VSTAR) is a modular code written in Fortran, which 
solves the line-by-line, radiative transfer problem for light in 
a stratified atmosphere with a defined chemical composition 
and a specified size distribution of cloud particles. The code 
is composed of several packages (MOD, LIN, RAY, PART 
and RT) which contain routines called by a user written 
program to build a specific atmospheric model. 




length (fiin) 



Figure 2. The observed spectrum (grey or red in the online 
version) of Jupiter's equatorial region (EZ) at K and H band 
overlapped with the model (black), which also cover the region 
between 1.85 and 2 (im characterized by the strong H2O absorp- 
tion in the Earth's atmosphere and thus it is difficult to observe 
from the Earth. The main discrepancy between the model and ob- 
servations appears at the low absorption region below 1.65 (im, 
which probes atmosphere levels below 1 bar. In addition we find 
typically larger mismatch in intensity of the methane absorption 
bands at H band than in K band, which is discussed in Section [4] 



The MOD package defines a two dimensional grid, 
which stores the values of parameters for the atmospheric 
layers in one dimension and the array of spectral points at 
which the model is calculated in the second dimension. The 
LIN package reads spectral lines from the line databases, cal- 
culate line profiles and sets the optical depth for each layer 
of the atmosphere. The RAY package calculates Rayleigh 
scattering optical depth for the atmospheric layers. The par- 
ticle scattering from a erosols and clouds is ca lculated with 
a Mie scattering code (jMishchenko et alJl20o3 l The PART 
package allows up to 20 different scattering particle modes 
with different size distributions, refractive indices and verti- 
cal distributions to be included in the model. The workhorse 
of the VSTAR code is the RT package, which performs the 
radiative transfer solution for each layer of the atmosphere 
using the data set up by the previously described packages. 
Different radiative transfer solvers can be selected but for 
the purpose of this wor k we use DISORT, the discrete ordi- 
nate solver described in IStanmes et all l|l98Sh . 

We model Jupiter's atmosphere using 36 vertical layers 
with specific physical properties given by pressure, temper- 
ature, as well as chemical composition expressed in terms of 
gas mixing ratios. The P-T profile of the planet's atmosphere 
in the pressure range between 0.001-1 bar has been stud- 
ied with infrared radiometry, radio occultations and direct 
sensi ng with the Galileo probe ASI experiment (ISeiff et ail 
1998). Most recent retrievals of atmosp heric temperatures 
from t he Cassini CIRS data published bv lSimon-Miller et al.l 
(2006) show time- variable, zonal differences in P-T profiles, 
which are also present in the Voyager IRIS data. In the range 
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Table 1. Jupiter observations on the 25 May 2007 



Band 


Centre 

H 


Resolving power 


Mean dispersion 
(nm/pixel) 


Start (UT) 
(h) 


Integration Time 

(s) 


Total 

(s) 


H 


1.637 


2270 


0.341 


14:40 


12 


48 


K 


2.249 


2250 


0.442 


14:58 


45 


180 



Table 2. List of selected regions of Jupiter 



Region 


1 


North Pole (NP) 


2 


Northern Upper Belt (NUB) 


3 


Northern Zone (NZ) 


4 


Northern Equatorial Belt (NEB) 


5 


Equatorial Zone (EZ) 


6 Equatorial Belt (EB) 


7 


Southern Equatorial Zone (SEZ) a 


8 Southern Belt (SB) a 


9 


Southern Zone (SZ) 6 


10 


Southern Upper Belt (SUB) 


11 


South Pole (SP) 



"These spectra are affected by a contribution of light reflected 
by Europa. b A spectrum of the atmospheric region with a bright 
'white' spot. 

of pressures probed with our model the strongest zonal vari- 
ations are of the order of a few Kelvin in the upper tro- 
poshere. 

In figure [4] we plot the av eraged profile from Voyager 
radio occultation experiments (|Linda]||l992t ). and the pro- 
file obtained from the Galileo probe descending into the 5- 
/xm h ot spot region in the north equatorial belt (|Seiff et al.l 
Il998l ). The notable differences in Voyager and Galileo pro- 
files appear above about 300 mb in the upper troposphere, 
which is expected since both profiles were probing different 
locations on Jupiter surface at different times. The 5-/xm 
hot spots are relat ively dry, localized regions with a low 
level of cloud cover (|Showman fc Dowlindl2000l ) . We believe 
the Voyager occultation data are likely to be more robust in 
describing the average properties of zones and belts consid- 
ered in our models. We tested the effects of changed tem- 
perature corresponding to the profile differences observed in 
Voyager ingress, egress occultations data and Galileo probe 
measurements for the Equatorial belt, Equatorial zone and 
the North polar region. The typical example of a difference 
between models with Voyager and Galileo derived P-T pro- 
files is shown in Figure [4] The differences are highest in the 
spectral ranges corresponding to the tropospheric pressures. 
However they are of the order of magnitude lower than the 
typical residuals of our fits to data, which are dominated by 
the opacities of different vertical regions in the atmosphere. 

Therefore for simplicity we used a single P-T profile to 
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Figure 3. The observed spectrum of Jupiter's equatorial region 
(EZ) at K and H band is shown with the difference between two 
models (in the lower panel) obtained with the Galileo Probe and 
the Voyager radio occultation P-T profiles, as presented in figure|4] 

model zonal composition of the atmos phere based on the 
Voyager rad io occultation observations (|Lindal et al.|[l98ll : 
lLindallll992l ). 

The chemistry of Jupiter's atmosphere is complex and 
involves a large number of atomic a nd molecular specie s pro- 
duced in photochemical processes (|Atreva et al. 1997). Our 
model considers only the major components of the atmo- 
sphere, which influence the shape of the observed spectra in 
the near infrared region. The most abundant atmospheric 
gases are H2 and He that produce smooth, collision-induced 
absorption features around 2.12 /jm. The other molecules, 
which absorb in the infrared spectral region, are methane 
(CH4) and ammonia (NH3). Methane is the most important 
absorber in the H and K band spectra. The mixing ratio of 
CH4 adopted in our model is 2.3 x 10 -3 , which is consis- 
tent with the latest in situ measur ements with the G alileo 
Probe Mass Spectrometer (GPMS) (IWong et alj|2004l ). The 
mixing ratio of met hane is relatively c onstant at pressures 
higher than 1 mbar ijMoses et al.ll2005h . 

3.1 The Methane Spectrum 

A tetrahedral molecule like methane has four vibrational 
modes, two bending modes described by quantum numbers 
v\ and v-j, and two stretching modes V2 and V4, so that 
each vibrational level can be described by its four vibra- 
tional quantum numbers (fi, 1% f3, va)- The frequencies of 
the two stretching modes are both about 1500 cm -1 and 
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Figure 4. Comparison between the temperature-pressure Voy- 
ager (solid line) and Galileo probe (dash line) profiles. 

the two bending modes both about 3000 cm -1 . These co- 
incidences result in a complex series of interacting vibra- 
tional states repeating at intervals of ~1500 cm -1 known 
as polyads. Polyads are generally designated by P n where 
n is the polyad number defined in terms of the vibrational 
quantum numbers by 

n = 2(i/i + u 3 ) + v 2 + vi (2) 

Po denotes the ground state. The first polyad Pi is also 
known as the Dyad since it contains the two vibrational 
levels (0, 1,0,0) and (0,0,0, 1). The P 2 polyad (or Pentad) 
contains 5 vibrational levels (1, 0, 0, 0), (0, 2, 0, 0), (0, 0, 1, 0), 
(0,0,0,2) and (0,1,0,1). Similarly P 3 is the Octad (8 lev- 
els) and P4 is the Tetradecad (14 levels). The number of 
levels grows rapidly for higher polyad numbers and results 
in complex systems of interacting bands which are increas- 
ingly difficult to model as the polyad number grows. 

A significant recent development is a new global anal- 
ysis of the methane spectrum over the wavenumber range 
0-4800 cm -1 (levels up to the Octad) reported by Albert et 
al. (2009). The new model fits line positions and intensities 
much better than prev ious analyses of individual polyads 
(e.g. iHilico et all l200ll ). The new model forms the basis 
of an improved set of methane line parame ters included in 
the 2 008 edition of the HITRAN database (|Rothman et all 
2009). It provides a good model of the low temperature 
methane spectrum for wavelengths greater than about 2.1 
/Urn. 

However, at shorter wavelengths (higher polyads) HI- 
TRAN 2008 lists mostly empirical line parameters from 
lBrownll|2005l ). These are measured line intensities and posi- 
tions at room temperature, and mostly lack quantum iden- 
tifications and lower-state energies. The lack of a lower-state 
energy means that line intensities cannot be reliably deter- 
mined for temperatures other than those of the original mea- 
surement. These lines are therefore not useful for modelling 



the spectra of the giant planets. However, some lines in the 
strongest part of the P4 band system have empirically deter- 
mined lower state energies derived from measuring the lines 
at tw o or more different tem pera tures. Thes e lines , with data 
from iMargolis et al.l (|l990l ) and lGao et~aH ({2009) cover the 
wavenumber range from 5500-6180 cm -1 . 

For this work, we therefore constructed a methane line 
list as follows. We used line data from HITRAN 2008 for all 
wavenumbers up to 4800 cm" 1 . We also used HITRAN 2008 
for the region from 5500-5550 cm" 1 , where empirical lower 
state energies are available. For the wavenumber range from 
4700-5500 cm -1 , we used line param eters calculated using 
the S pherical Top Data System (STDS lWenger fc Champion! 
1998). This range overlaps slightly with the range from 
HITRAN, but the lines from HITRAN are all in the P 3 
polyad system in this region, whereas the STDS lines are 
all from the P4 polyad, so there is no duplication. While 
the STDS effective hamiltonian model for the Tetradecad 
region is a preliminary one and does not give line positions 
that accurately match data, it does provide an adequate de- 
scription of the overall band shape, and provides the lower 
state energies that are missing from HITRAN in these re- 
gions . For the range 555 0-6236 cm" 1 we used line data 
from [Nikitin et all 1 20101 ) supplement e d wit h low temper- 
ature empirical data from IWang et al.l (|2010l ). We have not 
attempted to model wavelengths shorter than 1.6 /im. How- 
ever, new empirical methane line data for these wavelength 
region s a re starting to become availab le |Campargue et all 
l2010al lbl; ISciamma-O'Brien et alll2009l ). which may enable 
our modelling to be extended to shorter wavelengths. 

The line width data provided in HITRAN are for broad- 
ening in air. Data on broadening of meth ane lines in t he 1/3 
band in a variety of broadening gases bv lPind ||1992D , and 
iPine fc Gabardl feoosl ) shows that line widths broadened by 
H2 are about 2% higher than those broadened by N2, those 
in O2 are about 6% less, and those in He are about 35% 
lower. This suggests that line widths in an H2/He atmo- 
sphere should be very similar to those in an N2/O2 atmo- 
sphere. We have therefore used the air broadening data from 
HITRAN to determine the line widths in this work. Where 
the line data are taken from sources without broadening pa- 
rameters, we used the default line half width value of 7 = 
0.06 cm -1 atm" 1 and a tempe r ature exponent N = 0.85 as 
recommended bv lNikitin et al.l (|2010l ). 

3.2 Other Absorbers 

The atmospheric abundance of NH3 at high pressures 
was measured directly with GPMS and from the attenu- 
ation of radio signa l sent between the probe and orbiter 
(Folk ner et al.| [l998). The mixing ratio appears to decrease 
from 10~ 4 at pressures above 2 bar to 10~ 6 below 1 bar. 
Ammonia is known to exist in gaseous form and as con- 
densed clouds in Jupiter's troposphere up to pressures of 
0.1 bar l|Kunde et al.lll982D . but it is rapidly depleted in the 
stratosphere, due to UV photolysis. 

The mixing ratios of NH3 adopted in our model are 
constant (2 x 10~ 4 ) in layers between 1 and 0.5 bar. They 
decrease to 2 x 10" 7 at the 100 mbar along the saturated 
vapour pressure curve. The resulting profile i s consistent 
with m ixing ratios derived from the radio data in lSault et al.l 
|2004l ) above 0.5 bar. They note small differences between 
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the Equatorial Zone and North Equatorial Belt profiles, 
which we also tested in our models. As a result we observe 
changes in the absorption spectrum between 1.8 and 2 fim, 
where the v\ + va and vz + U4 bands of NH3 are present. We 
found that removing ammonia from the upper atmosphere 
(< 0.5 bar), which is most likely the case in down- welling 
regions, like belts and hot spots, does not have a significant 
effect on the K and H band spectra which are dominated by 
CH4 absorption. 

Absorption of ammonia was derived from the new list of 
Yurchenko et al. (2009) which contains over 3 million lines 
and is considerably more complete than the line parameters 
in HITRAN 2008, including improved shorter wavelength 
coverage. 

The bulk of Jupiter's atmosphere consists of H2 and 
He. Rotovibrational absorption bands due to collisions be- 
tween pairs of H2-H2 and ti2-He exist in the near infrared 
region under consideration, where they form a smooth fea- 
ture centred around 2.1^tm. Our model includes absorption 
coefficients for H2-H2 and H 2-He pairs at low temperatu res, 
derived by iBorvsowl (|2002i ) and iBorvsow et al.l (| 19891 ) re- 
spectively. 

It is worth noting that the 2.1/xm collision induced ab- 
sorption (CIA) depends also on the para-ortho ratio of hy- 
drogen spin isomers in the atmosphere of Jupiter. Our cur- 
rent models assume the equilibrium ratio at all tempera- 
tures. The Voyager IRIS data was used by IConrath et al.l 
( 1998) to demonstrate that the distribution of that ratio 
varies in the different zonal and vertical regions of the planet. 
The deviation from so called 'normal', 1:3, room tempera- 
ture, equilibrium value is relatively low for Jupiter, but it 
can produce a noticable effect on the shape of the 2.1/im 
feature in the spectrum. 

Our intention is to extend the capability of the VSTAR 
package in the near future to allow differences of the para- 
ortho H ratio as a function of atmospheric height, which will 
allow more accurate fit of the CIA absorption in atmospheres 
of cold giant planets, where the observed para-ortho ratios 
deviate even further from the normal value. 

3.3 Scattering 

In our model Rayleigh scattering from molecules and parti- 
cle scattering from aerosols suspended in the atmosphere 
was considered. The Rayleigh scattering at wavelengths 
above 1 /im is much reduced in comparison with scatter - 
ing from clouds. The Rayleigh cross sections jLioul 120021 1 
were calculated for scattering in a mixture of H2 and He, 
the dominant components of the atmosphere. 

The vertical structure of aerosols in Jupiter's atmo- 
sphere, derived directly from Galileo observat ions was sur- 
prisin gly different from theoretical predictions (|Atreva et al.l 
1 19971 ). which can be attributed to the unique region of the 
Probe entry. The Galileo Solid State Imager (SSI) observa- 
tions covered more extended region s over belts and zones 
as described in lBanfield et all (| 19981 ). The SSI images taken 
close to the limb of Jupiter gave good information about the 
clou d opacity in th e strato sphere. 

Banfield et al.l (|l998l ) used a retrieval technique 

ijBanfield et al.lll996l ) to obtain the scatterer density with 
altitude in the same spectral region as our observations. The 
upper troposphere between 700 and 100 mbar, and strato- 



sphere at pressures lower than 100 mbar defined by Banfield 
et al. are within the spectral sensitivity of the H and K 
bands. The results of the Banfield et al. retrieval of particle 
distribution from the Galileo imaging suggest two enhance- 
ments in cloud opacity corresponding to the upper haze be- 
tween 200 up to 20 mbar, and at the base of the tropospheric, 
condensation cloud at about 850 mbar. In equatorial regions, 
the upper haze is elevated to even lower pressures. 

The layer of smaller particles is believed to form a rel- 
atively uniform top of the tropospheric cloud, which con- 
tains larger particles of presumably NH3 ice and other con- 
densates. Observation s of limb darkening in the ultraviolet 
l|Tomasko et al.lll98fj ) constrain the size of the stratospheric 
haze to between 0.2 and 0.5 ^m. To model the clouds we 
chose a simple two-component distribution of mean parti- 
cle sizes, 0.3/im for the stratosphere and the conservative 
value of 1.3 fim size species for the upper tropospheric cloud. 
The combined observations from NIMS and NIR identify 
the lower tropospheric cloud as a mixture of two types of 
particles with 50 and 0.45 /im size, consistent with NH4SH 
molecule (|lrwinlll999l ). The upper troposphere, which our 
data is sensitive to, is believed to consists of 0.75 /jm size 
ammonia particles. 

The scattering properties of these clouds are calculated 
using Lorentz-Mie theory for a power-law particle size dis- 
tribution. The code calculates the asymmetry factor, extinc- 
tion and scattering coefficient, which are required in the 
radiative transfer calculations. The refractive index of the 
cloud particles is appr oximated by the r efract ive index of 
ammonia ice given by iMartonchik et all l|l984l ) for a wide 
range of wavelengths. 

In our model the vertical structure of the particle clouds 
is varied through the input list of optical depths for each 
layer of the atmosphere at the reference wavelength of 1.5 
^m. Finally the boundary conditions for the model are given. 
Solar illuminat ion at a given beam angle is provided us- 
ing data from iKurucj (|2009l ) and a non-reflecting surface 
is specified at the base of the atmosphere. The set of line 
data, particle properties and their distribution form inputs 
for the DISORT radiation transfer equation solver, which re- 
turns the fluxes and radiances at the top of the atmosphere 
for each wavelength specified in the grid. The models were 
calculated at a spectral resolution of 0.01 cm -1 , sufficient to 
resolve the spectral structure, and then convolved with an 
instrumental point spread function for comparison with the 
observed spectra. 



4 MODEL OF THE BANDS IN THE 
ATMOSPHERE OF JUPITER 

To show the functionality of the VSTAR code, we used the 
most recent observational data described in the previous sec- 
tion, which provided the input parameters for characterisa- 
tion of Jupiter's atmosphere. We tested the sensitivity of our 
models to changes in these parameters and confirmed that 
spectral features are dominated by methane absorption with 
the addition of clouds above the 800 mbar level. 

A simple model of the longitudinal belts and zones in 
the atmosphere of Jupiter was derived by manually vary- 
ing the amount and vertical distribution of cloud opacity at 
different heights, while the chemical composition of the at- 
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Figure 5. In the upper panels the observed spectra of the se- 
lected zones and belts in the K band (on the logarithmic inten- 
sity scale) are compared with the corresponding VSTAR models 
derived for these regions. In the lower panel the observed spec- 
tra in the H band are compared with the VSTAR models. The 
panels on the right (top and bottom) demonstrate the residuals 
obtained by taking the difference between the modelled spectra 
and corresponding data. 



mosphere was left constant across the surface of the planet. 
Since this is clearly not practical to explore full range of 
opacities at all levels of the atmosphere, we resorted to test- 
ing different combinations of opacities in the iterative way 
and fine-tuned them until a reasonable fit, with the smallest 
scatter in residual differences between spectra and models 
was achieved for both spectral bands H and K. 

Figure [5] shows the line- by-line models for the selected 



bands and zones with a resolution that matches the observed 
spectra. Our models mirror most spectral features remark- 
ably accurately considering that the average spectrum is a 
composite of sub-regions with possibly non-uniform cloud 
composition as evident in high resolution images from the 
Galileo SSI probes (Simon-Miller et al. 2001). We found that 
the most significant differences between the models and the 
data occurred in the equatorial zones, in spectral regions 
that correspond to the upper layers of the troposphere. 

The most recognizable features of the K band spectrum 
are the H2-H2 collision-induced broad absorption centered 
at 2.11 /im and stronger absorption in the 2.2 to 2.4 /im 
region due to methane bands. The narrower CH4 absorption 
features are due to the Q branches (A J = 0) of ^2 + ^3 at 2.2 
/im, + at 2.32 pirn and ^1+^4 at 2.37 fim. The latter two 
bands are not always clearly visible in the observed spectra 
due to low signal-to-noise in the high methane absorption 
region and a decreased sensitivity at the edge of K band 
filter. 

Spectra of the polar regions show the broad methane, 
and H2-H2 collision induced absorption much reduced as 
compared with the equatorial zones and belts. This can be 
understood by noting that the spectra of the polar regions 
sample relatively longer path-lengths through the higher lay- 
ers of the atmosphere, which unlike the deep layers con- 
tribute little to the H2-H2 absorption. 

The models are constrained by a simultaneous fit of ab- 
sorption levels for the broad and narrow band features of the 
spectra in both bands H and K. The lower panel of Figure[5] 
shows the spectra and models of the H band. The narrow 
spectral lines in the models tend to agree better with the 
observed spectrum than the broad shape and the intensity 
of absorption. We attribute this to the less accurate models 
of methane line parameters in this spectral region as dis- 
cussed in section |3T1 The accuracy of the wavelength fit is 
limited mostly by the spectral calibration of the data, dom- 
inated by the correction of spectral curvature in the IRIS2 
grism system. The wavelength of clearly identified lines in 
the data and models agrees within 5 x 10 _4 //m. Most dis- 
crepancies in models and observed spectra are due to the 
differences in intensity of a particular line or line series in 
methane absorption. 

The Jupiter K band spectrum is defined by the opaci- 
ties distributed between about 700 and 2 mbar. The shape 
of collision induced H2-H2 absorption, which samples deep 
layers of the atmosphere is highly sensitive to the cloud op- 
tical depth in the Jupiter's troposphere, while the level of 
broad methane absorption beyond 2.2 /jm varies dramati- 
cally in response to changes in the stratospheric cloud op- 
tical depth. The H band spectral region samples the higher 
levels of the troposphere. 

Table [3] shows the total optical depths of the tropo- 
spheric (1.3 /im particles) and stratospheric aerosols (0.3 
/im particles), which gave the best fit for each spectrum. 
All zones and belts show differences in the vertical distri- 
bution and the strength of opacity. Most notably the atmo- 
sphere appears to be relatively clear between 200 and 40 
mbar in the northern and equatorial regions, w hich is con- 
sistent with conclusions in lBanfield et al] (|l998l ). who find a 
similar clearing around 100 mbar in their models in the same 
spectral region. We find that no trial models with added 
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Table 3. Total cloud optical depths in zones and belts 



Particles 


1.3/im 




0.3/jm 




Zones 


Extent 
(mb) 


T 


Extent 
(mb) 


T 


NP 


563-631 


1.50 


56-200 

'3 7 


0.0100 

U.UUoU 


NUB 


316-562 


0.52 


125-200 
7-9 


0.0003 
n nni ^ 

U.UUIO 

0.0001 


NZ 


316-631 


1.72 


22-25 


0.0010 


NEB 


355-751 


5.70 


25-35 


0.0018 


EZ 


251-631 


5.72 


18-22 


0.0023 


EB 


316-708 


5.09 


25-35 


0.0070 


SEZ a 


447-708 


3.10 


39-50 


0.0720 


SB" 


355-708 


6.06 


141-251 
40-50 


0.0188 
0.0090 


sz b 


446-631 


2.25 


100-355 
25-35 
3-5 


0.0295 
0.0060 
0.0010 


SUB 


316-751 


1.29 


141-316 
25-40 


0.0024 
0.0080 


SP 


251-708 


2.20 


22-251 
2-4 


0.0038 
0.0012 



tems, which forms slightly higher clouds. Our model of the 
SZ white spot requires a relatively thick cloud at the 100- 
251 mb pressures with the small size particles, such as used 
for the stratospheric layers, which provided better fit to the 
observed spectra than any larger particles. 

The zones clear of Europa light show higher strato- 
spheric cloud opacities, which are distributed more thinly 
than clouds in nearby belts. The stratosphere in the high 
latitude northern (NUB) zone is more spread-out but rela- 
tively thin compared with the equatorial regions. The strato- 
sphere above the southern upper (SUB) belt is also vertically 
extended, but in contrast to the NUB, has higher opacity . 
Previous models l|West et al.lll986l : ISimon-Miller et aLll200lh 
suggested that downwelling of the haze in the belt regions 
could produce thinner cloud at high pressures. Our models 
show unambiguously increased opacity of high stratospheric 
layers in polar regions, which is consistent with the p revious 
studies of Jupiter's atmosphere ijPrvor fc Hordlfl99ll l. 

In Figure [6] the models of remaining regions are pre- 
sented. Northern and southern zones (NZ and SZ) display a 
slightly different shape of the collisional absorption, which 
suggests the existence of more compressed, optically thicker 
tropospheric cloud in SZ than the cloud in the NZ. The SEZ 
and SB spectra (Figure [TJl show the additional opacity due 
to Europa's spectrum dominated by fine-grained water ice 
(Carlson et al. 1996). The effects of this opacity are clearly 
visible in the spectral region beyond 2.15 /im, which is sen- 
sitive to the highest atmospheric levels in our model. 



"The spectra affected by a contribution of light reflected by 
Europa, which boosts reflection in the otherwise most highly 
absorbing methane bands regions in K and H bands. This 
additional reflection was not included in the modelling. 'A 
spectrum of the atmospheric region with a bright, 'white' spot, 
which can be modelled as the relatively thick stratospheric 
cloud extending upwards to 100 mbar. 

cloud opacity in this pressure range were able to reproduce 
the observed spectra. 

We find that the clouds tend to lie deeper in belts than 
in zones. The tropospheric cloud seems to extend up to 250 
mb in the equatorial region with the highest total opacity of 
the order of r = 5, while the stratospheric haze is very thin 
with the opacities below r ~ 0.005. Our models suggest that 
this haze may be slightly thicker in the southern hemisphere 
of the planet. However the SEZ and SB regions are affected 
by the spectrum of Europa, while the spectrum of the SZ 
zone may reflect a local cloud pattern. 

The image of Jupiter (Figure [TJ shows a rather large 
white spot in this region, which is subtended by a slit. The 
high resolution images taken by the New Horizons LEISA 
imager at 1.53 and 1.88 fim in February 200/Q 

show sim- 
ilar weather features in the region of SZ and SUB. This 
spot is most likely a storm system such as the anticyclonic 
white oval modelled in Simon-Miller et al. (2001). They note 
that fresh material is upwelled in such high-pressure sys- 

1 On-line images available from the NASA/ Johns Hopkins Uni- 
versity Applied Physics Laboratory/Southwest Research Institute 
\http://photojoumal.jpl.nasa. gov/catalog/PIA0925l\ 



5 DISCUSSION 

The line-by-line radiative transfer models for high resolution 
spectra from currently available ground-based instruments 
serve as an important test of the consistency in our un- 
derstanding of the composition and processes that govern 
planetary atmospheres. 

We achieved a good fit to the data assuming rela- 
tively simple physical structure and chemical composition 
of Jupiter's atmosphere that was kept uniform across the 
whole surface of the planet. The real atmosphere is cer- 
tainly more complex. Lacking a detailed P-T profile for each 
zone and belt we used a single profile of the atmosp here for 
all re gions based on radio occultation observations l|Lindall 
Il992h . This seems justified despite small differences in atmo- 
spheric profiles observed by Voyager 1 in frared spectrometer 
(IRIS ) and Cassini CIRS experiments (|Simon-Miller et al.l 
2006) for different locations on the planet. It is argued, that 
because horizontal differences in observed brightness tem- 
perature are small, the larger deviations in observed bright- 
ness te mperature are pro bably due to varied opacity in the 
clouds (jWest et al.lll986l ). 

We also used a relatively simple representation of the 
chemical composition of the atmosphere, despite the fact 
that the abundances of gases in different regions of the atmo- 
sphere may vary with height due to evidently complex, plan- 
etary weather system. The upwelling in zones may cause the 
depletion of volatiles due to condensation and precipitation. 
In the down-well regions the mixing ratios of such volatiles 
may be increased in deeper levels of the atmosphere. While 
methane is well mixed within the vertical region of our sensi- 
tivity, there have been zonal differences in mixing ratios pro- 
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Figure 6. The observed spectra, models and the difference be- 
tween them (on the linear intensity scale) are shown for the north- 
ern and southern, high latitude zones (NZ and SZ) and belts 
(NUB and SUB). Although the spectra of the northen regions 
show typically stronger absorption visible in both H2-H2 feature 
and the wavelength region above 2.2 jim, this differences can be 
explained in terms of 'local' weather effects as discussed in Sec- 
tion g] 



files o f ammonia below 600 mbar suggested by ISault et al.l 
l|2004h . Our models turned out to be not sensitive to these 
differences due to combined effect of the low NH3 absorp- 
tion in our spectral range and decreased sensitivity of our 
models to layers below ~ 700 mbar. We find that it is suffi- 
cient to consider only the major contributors to the spectral 
absorption at H and K bands: H2, He, CH4 and NH3 and 
to modify the opacities of aerosols at different heights. How- 
ever our systematically worse fits in equatorial zones may 
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Figure 7. The spectra of zones affected by opacity of Europa 
are shown compared with the simple models, which neglect the 
scattering contribution from this moon. The differences between 
spectra and models (bottom panels) show this clearly for spectral 
regions above 2.2/^m . 



reflect the unrecognized complexity in these regions, which 
is not described in sufficient detail by the initial parameters 
of our model, such as P-T profile and mixing ratios of the 
chemical components. 

The suggested constituents of clouds in the troposphere 
below 500 mbar are H2O, NH4SH and NH3. Ammonia con- 
densation occurs at the highest level at about 700 mbar at 
its saturation point. The actual condensation level can be 
lowered or lifted due to d ifferences in the molar abundances 
of NH^ IWesT et al . (1986) suggested that the cloud needs to 
have a mixed composition of particles with sizes between 3 
and 100 /im to explain high opacity observed at 45 /im. The 
larger particles are concentrated at the base of the cloud 
at pressures below our sensitivity limit. Consistently our 
models are most sensitive to changes in the particle sizes 
below 3 fim. We find that the tropospheric cloud in belt re- 
gions occurs at typically higher pressures compared to zones, 
which could be interpreted as the result of vigorous convec- 
tive downwelling and upwelling process. 

In contrast the stable haze of ~ 1/xm-size particles, mix- 
ture of ammonia ices and residual chromophores responsible 
for the coloration of the clouds exists above the lower, con- 
vective troposphere. In the warmer temperatures there is 
no condensation expected and consequently no cloud for- 
mation. The origin of the stratospheric haze is thought to 
be due to photochemical processes. 

This original pict u re of the layers above 700 mbar 
shown in West et al.l |l986) was developed further by 
Banfi eTd et al.l 1 1998! ). whose model includes a clearing of 
the haze at about 100 mbar pressures. They explain it by 
the aerosol coagulation at this level of stratosphere, which is 
the cause of the increased, gravitational fallout of heavy par- 
ticles. We observe such a clearing in models of most of the 
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northern zones and belts of the planet. Above the clearing 
we find the optically thin haze of reflective particles with a 
predominantly lower sizes, than particles in the troposphere. 
This haze extends to the low pressures of 2 mbar in the po- 
lar regions. The chemical composition of the stratospheric 
haze is still debated. It is also no clear if the polar haze ob- 
served at lower pressures has the same composition as the 
equatorial stratosphere. 



6 SUMMARY 

We have used the new, line-by-line, radiative transfer pack- 
age, VSTAR, to model the near infrared spectra of Jupiter. 
The models were produced for regions with different ab- 
sorption properties implied by varied cloud patterns in the 
atmosphere of the planet. With our improved database of 
methane absorption lines at low temperatures, we achieved 
a good fit to the spectra obtained with the AAT/IRIS2. 
We were able to reproduce conclusions suggested in previ- 
ous studies about the atmospheric conditions at pressures 
below 1 bar, which probe the upper troposphere and strato- 
sphere of the planet. We showed that ground-based spec- 
troscopy with large telescopes can significantly enhance the 
results achieved from space. With continued improvements 
in methane line parameters, we will be able to extend our 
model of the near IR spectrum of Jupiter and other giant 
planets and to cover full range of wavelengths between 1 and 
2.5 (jm with almost arbitrary resolution. This will make it 
possible to search for and measure weak spectral features 
due to trace constituents, that at present cannot be picked 
out of the forest of methane lines. These include, for ex- 
ample, absorption due to monodeuterated methane (CH3D) 
with lines at ~1.55 /jm that can be used to measure the 
D/H ratio of the giant planets with much higher accuracy 
than has been possible so far. 

The most appealing feature of the VSTAR package is its 
versatility. Although slower than correlated-k techniques, it 
allows spectral modelling of a wide range of objects, from the 
cool brown dwarfs to cold planets and their moons. These 
models can be calculated at any required spectral resolution 
and can therefore be used in conjunction with data from 
instruments now available on large ground-based telescopes 
that can provide near-IR spectra at resolving powers as high 
as 100,000. With the ever increasing speed of computing de- 
vices, it may be soon practical to use line-by-line techniques 
like VSTAR to achieve best fits to the high resolution spec- 
tra via automated retrival schemes. 

The methane line list used in this work (in HI- 
TRAN 2004 format) is available for download from 
http://www.phys.unsw.edu.au/~jbailey/ch4.html 
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